Layer-by-layer assembly of colloidal particles deposited onto the polymer-grafted 

elastic substrate 
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We demonstrate a novel route of spatially organizing the colloid arrangements on the polymer- 
grafted substrate by use of self-consistent field and density functional theories. We find that grafting 
of polymers onto a substrate can effectively control spatial dispersions of deposited colloids as a result 
of the balance between colloidal settling force and entropically elastic force of brushes, and colloids 
can form unexpected ordered structures on a grafting substrate. The depositing process of col- 
loidal particles onto the elastic "soft" substrate includes two steps: brush-mediated one-dimensional 
arrangement of colloidal crystals and controlled layer-by-layer growth driven entropically by non- 
adsorbing polymer solvent with increasing the particles. The result indicates a possibility for the 
production of highly ordered and defect- free structures by simply using the grafted substrate instead 
of periodically patterned templates, under appropriate selection of colloidal size, effective depositing 
potential, and brush coverage density. 
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Colloidal particles can self-assemble into a rich variety 
of highly ordered structures 1] on periodically patterned 
substrate |^], block copolymer scaffolds [3j], vesicle sur- 
faces with the opposite charge HL and at brush/air in- 
terfaces d, Q and liquid/liquid [3] or water/air Q in- 
terfaces. The major challenges in this field are how 
to assemble monodispersed colloids into highly ordered 
structures with well-controlled sizes and shapes, and 
how to achieve layer-by-layer growth of ordered struc- 
tures. Sedimentation is a simplest approach for col- 
loidal crystallization, however, usually leads to uniform 
or simple close-packed arrays of colloidal particles on 
smooth substrates. Effective control over interaction 
and arrangement of colloids [1, is possible by using 
densely polymer-grafted substrates since the entrop- 
ically elastic energy of brushes is comparable to thermal 
energy, and self-assembly depends critically on thermal 
energy j^]. A balance between depositing force of col- 
loids and entropic force of brushes probably leads to the 
formation of ordered colloidal crystal structures, rather 
than highly disordered aggregates. 

The grafting of polymers to surfaces is a simple and 
useful approach to stabilize colloids against aggregation 
and adsorption. It was experimentally reported [111 ] that 
the thickness of brushes can be adjusted between several 
nm and 1 ^m. The grafted polymer always exerts a repul- 
sive entropic force on incoming particles, and past works 
focused on the interaction between polymer brushes and 
individual incoming particles and how to prevent the ad- 
sorption of colloids such as proteins onto surfaces under 
various grafting density, chain length, and interactions 
between chains and surface ■ To the best of our knowl- 
edge, there have been no systematic theoretical studies 
into the self-assembling structures of colloidal particles 
when deposited onto the grafting substrate. Here, we 
undertake the first theoretical study of deposition of col- 
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loidal particles from polymer solvent onto the grafted 
substrate, and address an important issue about how to 
realize the spatial organizing dispersions of colloidal par- 
ticles. 

We consider a mixture of np solvent chains and Up 
colloidal particles with the radius R, in contact with a 
substrate grafted with tIq, polymer chains. The substrate 
is horizontally placed in a:y-plane which is positioned at 
2 = 0, and all polymer chains and colloids are allowed 
in the region Q < z < Zmax- Z^ax is large enough to 
avoid the influence of brushes and particles, namely the 
solvent chains attain their bulk property for large values 
of z It is shown that aggregates on the deformed 

polymer brush are anisotropic, and extend along only 
one direction for minimizing the entropy loss of brush 
chain conformations [6]. Here, the existence of "settling" 
force helps to form horizontally oriented cylinders. Thus 
we assume translational invariance along y for the sake 
of simplicity, and the calculation can be reduced to the 
xz plane. The volume of the system is x Zmax, 
where is the lateral length of the surfaces along the 
X axes. The grafting density is defined as a — Ua/Lx- 
All polymer chains are of the same polymerization in- 
dex N and flexible with the same statistical length a, 
and incompressible with a segment volume Pq^. The 
probability distribution for molecular conformations of 
a Gaussian chain a is assumed to take the Wiener form 
P[ra{s)] oc exp[-2f^ ds\^^^\^], where ra(s) denotes 
the position of segment s on chain a. Recently, Balazs 
and co-authors [1^ have successfully combined a self- 
consistent field (SCF) theory with a density functional 
theory (DFT) to study mixtures of diblock copolymer 
and nanoparticles. Here, we use the grand canonical form 
of SCF theory which has been proven to be powerful in 
calculatiiig equilibrium morphologies in polymeric sys- 
tem [13, [ij, [3, [III, to deal with polymer solvent and 
brushes, while particles are treated by DFT [l^, [l3| to 
determine their favorite distributions. The grand canon- 
ical partition function [l2j for the system can be written 
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where ^ is the chemical 



potential per solvent segment and Zn^^np,np is the canon- 
ical partition function for Ua grafted chains, n_g solvent 
polymers, and Up particles: 
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where the first 5 function enforces incompressibility and 
the second 5 function determines the position of the uni- 
formly anchored chain end. Ub is Boltzmann's constant, 
T is the temperature, and v is the interaction energy. 
Rp is the position of the center of the pth particle. The 
operators (ps, and (pp represent the local concentra- 
tions of grafted chains, solvent, and colloidal particles, 
respectively. The SCF theory gives the free energy F 
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where ip, (ps, and (pp are the local volume fractions of 
brushes, solvent chains, and particles, and the overall vol- 
ume fractions of brushes and particles are given by and 
Pp stands for the particle center distribution, and the 
local particle volume fraction is then given by </5p(r) — 
N- I\r'\<R'^^'Ppi^ + ^13]. Qc = /drg„(r,s)gj(r,s) 
represents the single chain partition function of brushes 
subject to the field w, and Qs = J drqs{r, s)qs{r,l — s) 
and Qp = / dr exp[— K;p(r)] are the partition functions 
for solvent and particles under fields Wg and Wp, respec- 
tively. The end-segment distribution functions qi{r,s) 
and g|(r, s) represent the probability of finding the s*'* 
segment at position r respectively from two distinct ends 
of chains, qi satisfies a modified diffusion equation ^ — 

^V^9i — 'Wi{r)qi, and qj meets the same diffusion equa- 
tion but with the right-hand side multiplied by —1. The 
last term in Eq.(2) is DFT term [iGj] accounting for the 
steric interaction between particles, and the excess free 
energy ^'(<p„) per particle is from the Carnahan-Starling 
function \l7j with the weighted particle density, ^p(r) 
p^ . The interaction energy v is given by Nv / p^ksT = 
J dr[xbsNip(ps + XbpNipipp + XspNipsifip + geZipp], where 
the x's are the Flory-Huggins interaction parameters be- 
tween the different chemical species. We fix N — 100, 
XbsN = 0, and XbpN — XspN = 12.0, since we assume 
that all the polymers have the same chemical nature, 
while the particles are insoluble to polymers [l^. In ad- 
dition, a depositing force is applied normal to the sub- 
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FIG. 1: Particles density distributions in x-z cross-sections 
under different particle volume fractions, (a) (jjp — 0.03, (b) 
(t>p = 0.04, (c) = 0.05, (d) (t)p = 0.06, (e) (/>p = 0.092, (f) 
(l)p = 0.1, and (g) (j>p = 0.145. The color scale bar shows the 
local density values of particles in Figs. l(a)-l(g). 



strate, and > is the strength of settling field act- 
ing on particles. Here, we choose ge within the range 
0.1 ^ 0.6 which, on one hand, can ensure that the size 
of particles(2i? — 0.5i?o) is comparable to the sedimen- 
tation length (oc 0.8i?o b.ORo, where 
Rq = aN^/'^ characterizes the natural size of polymer. 



And on the other hand, the competition between deposit- 
ing energy of colloids and elastic entropy of deformed 
brushes can significantly appear within this range of g^- 
The volume fraction (j)p of particles changes from 0.03 to 
0.15 for ensuring the possibility of colloidal crystallization 
on the brush surface and layer- by-layer growth of colloids. 
For sub-micrometer particles, the depositing force may 
be introduced either by gravity or by centrifuge [H, |20| . 
whereas for nano-sized colloids, the sedimentation may 
be controlled via electrophoretic techniques [13, HH, H^. 
We take cr — 0.25 so that particles hardly embed into the 
dense brush. A periodical boundary condition for the x- 
direction is applied, while for the z-direction, the region 
of z < is forbidden and the region z > Z^^ax is treated 
as a bath of solvent. Zmax is fixed to be 60a, and the lat- 
eral length Lx is selected to minimize the free energy of 
system [23j. All the sizes are in units of a. In SCF theory, 
the fields and densities are determined by locating saddle 
points in the free energy subject to the incompressibility: 
<p(r) -I- (/Js(r) + '/'p(r) = 1. The resulting SCF equations 
are solved by the combinatorial screening algorithm of 
Drolet and Fredrickson [T5| . 

We first examine the effects of colloid volume fraction 
(j)p on colloidal dispersions at the top of brushes. Figure 1 
shows the morphology changes of dispersed particles with 
increasing (j)p for the parameter ge = 0.5. We find that 
at low particle volume fraction (Fig. 1(a)), the dispersion 
of particles forms a laterally uniform layer on the top of 
brushes. This means that the deposition of particles does 
not deform brushes to in turn react with the dispersion of 
particles, and particles only spread over the grafting sur- 
face, which is similar to that of particles depositing onto 



FIG. 2: (a) The minimum spacing dc between cylinders vs 
Qs with error bars, (b) The critical correlation factor Ac and 
penetration depth Sc as a function of ge- 



the smooth substrate. As more colloidal particles are 
deposited, the grafted polymer is compressed under de- 
positing potential, and it responds with a restoring force 
which further drives particles to self-assemble into cer- 
tain structures for counteracting deposition of particles, 
in the requirement of minimizing combinational contri- 
bution of colloidal depositing energy and brush entropy. 
Therefore, colloidal particles assemble into colloidal crys- 
tals under sedimentation, and morphology of regularly 
separated cylinder structures emerges. The number of 
cylinders will increase with the volume fraction (f)p (Figs. 
l(b)-l(d)). However, further increase of (j)p leads to the 
formation of second layer of cylinders(Figs. 1(e) and 1(f)) 
piled on the first layer created by brushes, and even forms 
the third layer of cylinder structures (Fig. 1(g)). 

The first-layered colloidal structure is formed, due 
to the competition between colloidal sedimentation and 
brush entropy effects. For a fixed the deposited num- 
ber of cylinders of radius R increases with increasing (j)p 
when the variation of depositing potential per cylinder 
geSTrR^kBT /N > AFf,. Here, AFb represents a brush 
stretching energy penalty of a local single cylinder with 
the penetration depth e into the brush [2j , which is given 
by AFb = ^Tje^TT f25'|. For melt brush, the shear mod- 
ulus 1] — i]o — Sa^ksT. We introduce a modified fac- 
tor A (i.e., T] = A770) to account for the correlation effect 
due to the interplay between multi-cylinders mediated by 
brushes. Therefore, the variation of A will mainly depend 
on the distance d between cylinders. However, we should 
point out that here, the density profiles of free ends of 
brushes are diffusive due to the existence of solvent chains 
above brushes. This will lead to the great decrease of 
shear modulus compared to that of melt brushes. We 
have A < 3^"^^ , and the equality is valid at the criti- 
cal (j)p where the minimum spacing between cylinders is 
reached, and larger (/)p will lead to second-layer aggre- 
gation of particles. From the SCFT/DFT calculations, 
we determine the minimum spacing dc (Fig. 2(a)) and 
the critical penetration depth Sc of cyHnders(Fig. 2(b)), 
and find that dc decreases, whereas £c increases with ge- 
This clearly shows that with increasing g^, depositing po- 
tential of particles is balanced by further deformation of 
brushes. Figure 2(b) also gives the critical modified fac- 
tor Ac as a function of ge, showing that at critical (f)p, 



FIG. 3: (a) Lateral statistical concentration morphology 
of brush(dotted curve), free chains(solid curve), and parti- 
cles(dashed curve) for (ftp — 0.1. (b) Entropies of brushes (St) 
and solvent (Ss) vs 0p. 



the effective shear modulus first increases with increas- 
ing ge and then slightly decreases at large ge- For small 
ge, the shear modulus increases with ge, indicating that 
the correlation between cylinders is enhanced due to the 
decrease of d^. In contrast, for large ge, the minimum 
spacing between cylinders keeps almost unchanged, but 
the embedding depth Sc becomes larger than the cylinder 
radius so that the grafted chains can cross the narrow 
gap between cylinders and fill the upper space. Thus, 
the grafted chains are slightly released, which accounts 
for small decrease of Ac at large ge{> 0.5). 

The formed first-layer colloidal crystal may serve as 
a template for next-layer structural formation, and thus 
provides a possible route to the fabrication of multi-layer 
microstructures. We find from Fig. 1(f) that the second 
layer is well arranged, based on the already deposited 
layer. However, the morphology selection of second-layer 
colloidal dispersions will be affected by the entropic effect 
of solvent polymers. By calculating the z-direction aver- 
aged density profiles of brushes, free chains, and particles 
for the case of Fig. 1(f), Fig. 3(a) clearly shows that more 
solvent chains will fill the space between the first layer 
and the second layer of cylinders. Therefore, the second- 
layer cylinder structure is selected to alternately arrange 
with the first layer for increasing the configurational en- 
tropy of confined solvent chains. In fact, the alternating 
arrangement of cylinders in Fig. 1(g) further supports 
our viewpoint on controlled layer-by-layer growth driven 
entropically by polymer solvent. Figure 3(b) gives the 
entropies of brush and solvent as a function of (j>p. We 
find that when (t>p takes the range 0.09 ~ 0.10 corre- 
sponding to the forming process of second-layer struc- 
ture, the brush entropy retains almost unchanged, while 
the solvent entropy sharply declines, meaning that the 
second-layer colloidal assembly is out of brush effects, in- 
stead the solvent entropy dominates the final equilibrium 
dispersion of the second-layer particles. 

Finally, Fig. 4 shows the entropy of brushes and 
the root-mean-square fluctuation A of statistical brush 
height h as a function of ge for (pp — 0.06 and 0.1. 
We see that the brush entropy decreases with an increase 
of ge , but the height fluctuation due to brush deformation 
increases with ge- Correspondingly, the particle distribu- 
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FIG. 4: The brush entropy St and the modulated fluctuation 
A of brush height vs Qe- {a.)(l>p = 0.06. (b) (j)p = 0.1. 

tions in the inset of Fig. 4 signify the formation of cyhn- 
der structures with varying g^. Figure 4(a) shows that 
for small g^, the deposition of particles did not deform 
brushes which therefore do not react with the dispersion 
of particles. On the other hand, the relatively large range 
of the parameter g,, can stabilize the one-layer cylinder 
structures due to strongly entropic restoring forces of 
brushes. In contrast, Fig. 4(b) shows that a small range 
of ge may retain the two-layer cylinder structures. When 
ge is small, there are not enough particles deposited onto 
the top of brushes, leading to one-layer structure. As ge 
is relatively large, the two-layer structure is destroyed, in- 
stead the alternating structure of one- and bi-layer cylin- 
ders appears, because the entropy restoring force of sol- 
vent chains is weaker than that of brushes, and may not 
completely offset the settling energy of particles if the 



second layer is formed. It is actually interesting that the 
brushes have large entropic restoring forces which easily 
stabilize colloidal dispersions. For example, depending 
on the colloidal weight and volume fraction, brushes can 
adjust the number of cylinders formed, in contrast to col- 
loidal crystallization in non- adsorbing polymer solvent. 

In summary, we have demonstrated that under suitable 
density and depositing force of particles, colloidal par- 
ticles can be sorted into alternating arrays of cylinders 
by use of grafted substrates. The colloidal dispersions 
are dominated by the requirement of minimizing com- 
binational contribution of depositing potential of parti- 
cles and entropic restoring force of the deformed brushes. 
With an increase of colloidal additions, controlled layer- 
by-layer growth is driven by entropic effects of solvent 
chains. The advantage of the present approach is that 
control over arrangement of colloids did not rely on other 
patterned and phase-separated copolymer Q tem- 
plates but was achievable via a polymer-grafted substrate 
which is easily manufactured. The approach that under 
sedimentation, polymer entropic restoring force drives or- 
dered structure formation, will offer a simple and power- 
ful alternative for producing 2D and even 3D structures, 
and may open up an unexplored route for engineering 
highly ordered structures from colloidal building blocks. 

This work was supported by the National Natural Sci- 
ence Foundation of China under Grant Nos. 10334020, 
10021001, and 20490220. 



[1] G. W. Whitesides and M. Boncheva, Proc. Natl. Acad. 
Set. 99, 4769 (2002); J. Zhu et al.. Nature 387, 
883(1997). 

[2] A. V. Blaaderen, R. Ruel, and P. Wiltzius, Nature 385, 

321(1997); K. Lin et aL, Phy.s. Rev. Lett. 85, 1770(2000); 

Y. Yin et al., J. Am. Chem. Soc. 123, 8718(2001). 
[3] W. A. Lopes and H. M. Jaeger, Nature 414, 735(2001); 

M. J. Misner et al.. Adv. Mater. 15, 221(2003); M. R. 

Bockstaller et al., J. Am. Chem. Soc. 125, 5276 (2003). 
[4] L. Ramos et al, Science 286, 2325(1999). 
[5] Z. Liu et al.. Nana. Lett. 2, 219(2002). 
[6] J. U. Kim and B. O'Shaughnessy, Phys. Rev. Lett. 89, 

238301 (2002). 
[7] Y. Lin et al.. Science 299, 226(2003). 
[8] H. Wickman and J.N. Korley, Nature 393, 445(1998). 
[9] K. Chen and Y. Ma, J. Phys. Chem. B 109, 17617(2005); 

C. Ren and Y. Ma, J. Am. Chem. Soc. 128, 2733(2006). 
[10] S. T. Milner, Science 251, 905(1991); J. Satulovsky, M. 

A. Carignano, and I. Szleifer, Proc. Natl. Acad. Sci. 97, 

9037 (2000). 

[11] M. Biesalski and J. Ruhe, Macromolecules 32, 
2309(1999); J. Habicht, M. Schmidt, J. Ruhe, and D. 
Johannsmann, Langmuir 15, 2460(1999). 

[12] P. G. Ferreira, A. Ajdari and L. Leibler, Macromolecules 
31, 3994 (1998). 

[13] R. B. Thompson et al, Science 292, 2469 (2001); Macro- 
molecules 35, 1060 (2002). 



[14] M. W. Matsen and M. Schick, Phys. Rev. Lett. 72, 
2660(1994); D. Petera and M. Muthukumar, J. Chem. 
Phys. 109, 5101(1998); F. Schmid, J. Phys.:Condens. 
Matter 10, 8105 (1998); M. W. Matsen and F. S. Bates, 
J. Chem. Phys. 106, 2436(1997); P. Maniadis et al., 
Phys. Rev. E 69, 031801(2004); M. MuUer, Phys. Rev. 
E 65, 030802(2002); T. Geisinger, M. MuUer, and K. 
Binder, J. Chem. Phys. Ill, 5241(1999). 

[15] F. Drolet and G. H. Fredrickson, Phys. Rev. Lett. 83, 
4317 (1999); Macromolecules 34, 5317 (2001). 

[16] P. Tarazona, Mol. Phys. 52, 81(1984). 

[17] N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 
635(1969). 

[18] The interaction of colloidal particles can be adjusted by 

chemically coating ligands onto the surfaces of them. 
[19] T. Biben, J. P. Hansen, and J. L. Barrat, J. Chem. Phys. 

98, 7330(1993); C. I. Addison, J. P. Hansen, and A. A. 

Louis, ChemPhysChem. 6, 1760(2005). 
[20] M. Holgado et al., Langmuir 15, 4701(1999). 
[21] M. Trau, D. A. Saville, and LA. Aksay, Science 272, 

706(1996); M. Giersig and P. Mulvaney, J. Phys. Chem. 

97, 6334(1993); Langmuir 9, 3408(1993). 
[22] N. V. Dziomkina and G. J. Vancso, Soft Matter 1, 

265(2005). 

[23] Y. Bohbot-Raviv and Z. G. Wang, Phys. Rev. Lett. 85, 
3428(2000). 

[24] We define e = max(/i(a:)) — min{h{x)), where the local 



5 



brush height h{x) = 2J2^ 'fii^, z)z/ z)). 
[25] G. H. Fredrickson et al., Macromolecules 25, 2882(1992). 



